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Abstract — We study delay tolerant networking (DTN) and in 
particular, its capacity to store, carry and forward messages so 
that the messages eventually reach their final destination(s). We 
approach this broad question in the framework of percolation 
theory. To this end, we assume an elementary mobility model, 
where nodes arrive to an infinite plane according to a Poisson 
point process, move a certain distance £, and then depart. 
In this setting, we characterize the mean density of nodes 
required to support DTN style networking. In particular, under 
the given assumptions, we show that DTN communication is 
feasible when the mean node degree v is greater than 4 • r] c (-f), 
where parameter 7= £/d is the ratio of the distance I to the 
transmission range d, and tijci) is the critical reduced number 
density of tilted cylinders in a directed continuum percolation 
model. By means of Monte Carlo simulations, we give numerical 
values for r\ c (7). The asymptotic behavior of r\ c (7) when 7 tends 
to 00 is also derived from a fluid flow analysis. 

Index Terms — DTN, capacity, percolation, criticality, mobility 

I. Introduction 

Delay-tolerant networking (DTN) is a last-resort networking 
paradigm for mobile nodes when direct and multi-hop con- 
nections are infeasible, i.e., the network is not connected and 
nodes have to carry the messages to the next node (store-carry- 
forward). We study the capacity of DTN to deliver a message 
to its destination(s) in the framework of the percolation theory. 
To this end, we assume an elementary mobility model where 
nodes arrive according to a Poisson point process on a plane, 
move a certain distance i, and then depart. Nodes perform 
epidemic routing, i.e., whenever two nodes meet, all messages 
are exchanged. We are interested in finding the circumstances 
under which the lifetime of a message becomes infinite (with 
some positive probability) so that the message could reach a 
recipient located at an arbitrary distance from the source. We 
obtain a fundamental criticality condition that characterizes the 
sufficient mean density of nodes to this end. The criticality 
condition takes form v > 1^(7) = 4 ^(7), where v denotes 
the mean node degree (number of neighbors), 7 is the ratio 
of I to the transmission range d, and rj c (^) is an unknown 
function, which we determine in this paper. 

For a mobile DTN, we find that v c (j) < 1.52 for all 7, 
while only at v ss 4.51 a (non-DTN) wireless ad-hoc network 
percolates and a gigantic connected cluster emerges (T). In 
other words, DTN communication is possible over a very 
sparse network provided that the network's topology changes 
in time. A convenient characterization in fact is to say that 
DTN is a network that is super-critical in space-time. In 



practice, DTN is often sub-critical at a random time instant 
and the messages find their way to the destination through the 
store-carry-forward routing in space-time. 

This work is motivated by different opportunistic network- 
ing schemes. One such scheme is Floating Content, where 
nodes replicate messages only within the area where each 
message is deemed relevant. A criticality condition, charac- 
terizing the circumstances under which the lifetime of the 
message can be expected to be long, was established in |2] 
under the assumptions of mobile nodes and point contacts (i.e., 
the transmission range is small compared to the dimension 
of the area). Similarly, Beachnet |3j the aim is to produce 
quasi-periodic "information waves" carried by an underlying 
field of immobile nodes (say, devices on a beach). Beachnet 
is designed to offer regular content dissemination when nodes 
are (mostly) stationary and their collective cooperation, based 
on simple local transmission rules, is sought. Similar con- 
cepts for one-to-many content dissemination include hovering 
information [4] and ad-hoc podcasting (5), but our model 
naturally covers the special case of one-to-one messaging. 
These types of DTN schemes introduce many interesting 
problems. Perhaps the most fundamental question is whether a 
network can transport messages to their intended destinations 
or not, and this is also our focus in this paper. 

A. Background and related work 

We assume epidemic routing which was proposed by Vahdat 
and Becker in [6|. The operational principle is very simple: 
whenever two nodes meet they exchange the messages only 
one of them has. The performance analysis of the epidemic 
routing often assumes exponentially distributed (i.i.d.) inter- 
meeting times leading to Markovian models [7|. As the size of 
the network grows, solving Markov chains becomes infeasible, 
and Zhang et al., in |8 1 obtained ODEs as a fluid limit of such 
Markovian models. In above work, the spatial dimension has 
been abstracted away. 

Jacquet et al. [9] consider the propagation speed of the 
information in DTN for a fixed set of nodes moving in a 
finite region. The spatial dimension is explicitly present in 
their formulation. Somewhat related, Grossglauser and Tse 
(lOj studied the performance of adhoc networks under mo- 
bility. They focused on connected networks, where nodes can 
either communicate directly or via multi-hop connections, and 
showed that mobility increases the overall capacity in the 
network given the users tolerate additional delays. 



2 



In this paper, we consider a fundamental feasibility problem 
of DTN-style communication with the elementary mobility 
model by means of percolation theory fTT) , fl2) . Percolation 
theory has been succesfully applied to study the performance 
of wireless multi-hop networks since the early work by Gilbert 
|jT3) . Similarly as in JT3| , most of the work utilizes undirected 
planar continuum percolation models to argue, e.g., about the 
network's connectivity or capacity. 

In contrast, in fl4) , we assumed stationary nodes and 
showed that opportunistic content dissemination schemes, such 
as the floating content, can be analyzed by using a three- 
dimensional continuum percolation model. In this paper, we 
adapt the same approach, but instead of approximating the 
process by undirected percolation model, we consider the 
actual directed percolation characterizing the dissemination of 
messages by mobile nodes exactly in a DTN network. There 
are less results available for the directed percolation models 
than there are for the normal undirected percolation, and 
results for directed continuum models are even more scarce. 
The basic directed cases are the bond and site percolation 
models, where, e.g., bonds have fixed directions (e.g., up and 
right in a square lattice) [11], [ 15 1. Directed models have been 
also studied in the context of scale-free networks fl6| , which 
arise in Internet, social networks etc. 

The direction property can be defined many ways. The 
model considered in this paper, is a specific case of directed 
continuum percolation, where three-dimensional objects may 
belong to a given cluster fully or partially. In particular, the 
third dimension is time and, due to the causality, only the 
post-contact part of an object belongs to the given cluster. 
In this setting, we obtain an elegant fundamental result for 
the minimum feasible mean node degree to support DTN 
communication as a function of the ratio of movement to 
the transmission range. This critical percolation threshold is 
determined by Monte Carlo simulations. The obtained results 
define the feasible operation regime for large DTN networks. 

The rest of the paper is organized as follows. In Section |H| 
we introduce the model and the notation. Section Hill contains 
the analysis and the description of the Monte Carlo simu- 
lations. The numerical results are given in Section IV and 
Section [V] concludes the paper. 

II. Model and Notations 

In this section, we describe our model, including the arrival 
process, nodes' mobility and the assumptions regarding the 
radio communication. Throughout this paper, we implicitly 
assume a large network and a long distance, either in space 
and/or time, between the source and the destination. 

1) Node mobility: We assume that nodes arrive according 
to a Poisson point process on an infinite plane with rate density 
denoted by A [node/m 2 /s]. Mobility is specified with a constant 
movement I, i.e., during the constant time t a node moves 
the distance of £ to a random direction (isotropic). After the 
movement, the node departs from the system. Note that this 
activity cycle from "an arrival" to "a departure" can model an 
activity period of a node (from active state to inactive state), 
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Figure 1: Node mobility in space-time: (left) during time t, a node 
moves a distance of £ — L to a random direction in (x, y)-plane and 
departs. ( right) The volume covered corresponds to a tilted cylinder. 



or the departure can correspond to an event where a node 
replaces the information content with another. 
The node density in plane is 

h = A • t, 

The arrival-departure cycle can be interpreted as the time when 
a mobile node is participating in the DTN activity. In the 
basic stationary case, I = and the nodes remain still for 
the constant time t before departing. We study the capacity of 
a large network to transport messages in space-time. Thus, the 
node mobility model applies to the intermediate nodes, while 
the source and destination node can be thought to be, e.g., two 
permanent nodes located far apart from each other. 

2) Information dissemination: Whenever two nodes are 
within each others' transmission range d, i.e., the distance 
between them is less than d, they immediately exchange the 
messages. The assumption of a fixed transmission range is 
often referred to as the Gilbert's disc model or boolean model 
p3| . The mean number of neighbors a node has, i.e., the 
mean node degree is thus v — h ■ nd 2 . In the percolation 
theory, two objects are connected if they overlap and thus the 
corresponding radius r is half of the transmission range, 

r = d/2. 

3) Space-time: The movement and the radio coverage of a 
node in space-time, comprising (x, y)-plane and the time axis, 
are depicted in Fig. [T] The radio coverage at any time instant 
is a disc. As nodes move in time, the coverage in 3D space- 
time is a sheared cylinder. The tilt of the cylinder is defined 
by the ratio of the movement to the transmission range, 

a £ = L 
7 d 2r' 

Somewhat counter-intuitively, the number density of nodes 
(cylinders) in space-time is equal to the arrival rate A[j] 

n = X. 

Instead of number density, the density of objects is often 
expressed using either the reduced number density rj, 

rj = n ■ V, 

where V is the volume of the shape, or the volume fraction 
4>, for which it holds that <f> = 1 — e~ v . 

'Consider a volume V = AT, with area A in (x, y)-plane and height T. 
The mean number of nodes in V is XV , i.e., the node density n is A. 
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In our case, V — irr 2 ■ t independently of 7. Without lack of 
generality, we can scale the space-time. The time axis can be 
scaled by constant a, e.g., so that t = 1. Similarly, the (x,y)- 
plane can be scaled by constant (3, e.g., so that 2 • £ + d = 1. 
The ratio 7 is clearly invariant under such scalings. Similarly, 

n (x a - x fi- 2 and V oc a/3 2 , 

and thus also the reduced number density r\ = n ■ V and the 
volume fraction (f> are invariant. 

III. Analysis 

In this section, we first study the corresponding directed 
continuum percolation model characterizing the evolution of 
a message in space-time. Then we describe Monte Carlo 
simulations to estimate the critical percolation threshold, and 
also derive an asymptotic scaling law for the threshold. 

A. Directed continuum percolation 

We consider a dynamic system in space-time, where each 
node corresponds to a (tilted) cylinder as depicted in Fig. [T] 
Two nodes can communicate from the moment the distance 
between them becomes less than the constant transmission 
range d. It follows that the process describing the evolution 
of a message in the infinite plane is equivalent to a directed 
continuum percolation of aligned or sheared cylinders in three 
dimensions fl4) . The radius r of the cylinders with respect to 
percolation is equal to d/2, and the constant holding time t 
corresponds to the height of the cylinders. The moment two 
r-cylinders touch each other, messages can be transmitted. 
The causality imposes that the information flows only in the 
direction of the positive time axis. 

Formally, let a 6 denote that a path from node a to 
node b exists in space-time in the sense that a could send a 
message to b (potentially via a multi-hop connection). Due to 
the causality, is not a symmetric relation, 

a — * b =£> b — * a. 

Moreover, is not even transitive, 

a b and b c =£> a -*+ c, 

as a message from a to b may arrive too late for b to deliver 
it to c. Let S(a) denote the set of nodes reachable from a, 

S(a) = {b : a ^ b}. 

We are interested in the size (i.e., the cardinality) of S(a), and 
let random variable S, cluster size, denote this quantity, 

S=\S(a)\. 

To illustrate the directed percolation, consider a sample 
clusters depicted in Fig. [2] obtained without mobility (£ = 0). 
The time axis points towards the top-right corner and without 
mobility all cylinders are aligned with the time axis. The 
dark red node (cylinder) indicates the source of a message, 
light yellow color indicates that a node has not obtained the 
message yet, and green color means that the node has it (at 
the given time). Those nodes that never obtained the message 
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Figure 2: Sample realizations of dissemination of a message with 
stationary nodes. The reduced number density is r\ = 0.37. The red 
cylinder is the source, the light yellow color indicates that a node is 
waiting for the message and the green color that a node has acquired 
it. Nodes which did not acquire the message have been omitted. 




Figure 3: Sample realization with tilt 7 = 1 and 77 = 0.14 
comprising 65 nodes. 

have been omitted for clarity. Fig. [3] illustrates the situation 
with mobile nodes. From the figures we observe how the 
information "flows" only in the positive direction of time. 

When the number density of the cylinders (in space-time) 
is above the critical threshold, denoted by n c , there is a 
positive probability that a random cylinder belongs to an 
infinite cluster. Mathematically the critical percolation density 
is defined as, 

n c = inf{n I P{S = 00} > 0}. 

From the information dissemination point of view, this means 
that there is a positive chance that a given message gets 
distributed "everywhere" and does not go into extinction. 

The critical number density n c obviously depends on the 
dimensions of the cylinders and the amount of the movement t. 
Therefore, the critical percolation threshold is often expressed 
using either the critical reduced number density r\ c or the 
critical volume fraction </> c , which are both invariant to scaling 
as explained earlier. That is, for the aligned cylinders with 
7 = 0, it holds that r\ c and cb c are independent of the 
dimensions of the cylinder fl4) . In our case, the cylinders are 
tilted to a random direction, and consequently, the criticality 
threshold depends on the tilt ratio 7, and we have r\ c — r/ c (jy) 
(i.e., 7 is a shape parameter). Formally, 

r? c (7) 4 inf{7? I P{5 7 = ex)} > 0}. (1) 

Two larger sample clusters are shown in Fig. [4] The cluster 
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Figure 4: Sample realizations of a message dissemination with stationary nodes. Left figure illustrates a finite cluster with 969 nodes 
obtained with i] = 0.37. In the right figure, the number density is higher, r\ — 0.39, which is above the percolation threshold and the sample 
realization (presumably) percolates. 



in left is finite and the message goes to extinction after some 
time, while the cluster on the right appears to lead to an infinite 
cluster, i.e., the system percolates implying r\ > rj c . 

The volume of the (tilted) cylinders is V = irr 2 ■ t and the 
number density in space-time is n = A. Thus, r] = irr 2 ■ Xt, 
yielding 

z/ = 4ry. (2) 

Consequently, in terms of the mean node degree, the criticality 
condition states 

Ml) = 4t7 c ( 7 ). (3) 

An infinite cluster with respect to the directed percolation is 
naturally an infinite cluster also with respect to the undirected 
percolation. In [14], we already determined the critical (undi- 
rected) continuum percolation threshold for aligned cylinders, 
f]* = 0.3312(1), corresponding to the case £ = 7 = 0, which 
thus serves as a strict lower bound for the directed percolation 
in the same setting, rj c {0) > r/*(0) = 0.3312(1). 

When a DTN network percolates in space-time, it means 
that there is a positive probability that a message makes its 
way to the destination no matter how far the destination 
is (in space). In this sense, the criticality condition |3]l is 
a fundamental result for the transport capacity of DTN. At 
the same time, it gives no guarantees on the probability that 
this happens for a particular message (cf. strength of the 
percolation), or on how long it takes (which obviously depends 
on the distance among other things, cf. (9)). 

B. Methodology 

An elegant way to determine the percolation threshold is 
based on the asymptotic behavior of the cluster size S. In 

2 What happens when r\ = rj c is not known in general, but for r] < r] c , by 
definition, no infinite cluster emerges. 



particular, the tail behaves according to 

¥{S> S \v}~As 2 - T f{{r 1 -r lc )s a ), 

where r and a denote the so-called universal exponents and 
A is some (non-universal) constant. In three dimensions fT7) , 

t = 2.18906 ±0.00006, 
a = 0.4522 ± 0.0008. 



(4) 



Near the percolation threshold, where 77 w r/ c . The Taylor 
series for f(x) is f(x) = 1 + Bx + . . ., which gives 

P{S>s\r,}- s T - 2 ~ A + AB(rj - ri c )s a + ..., 

i.e., the quantity on the left-hand side becomes a constant when 
i] = rjc- The critical r/ c is then determined as follows |14|, 



1 18]-[20|. First one collects a large sample set of cluster sizes 
and stores the results to bins so that the fcth bin, denoted by 
Bk, corresponds to the number of samples with cluster size 
s > 2 k . For rj = rj c , the quantity B k ■ (2 k ) T ~ 2 then has a 
constant tail. Thus, by trial and error such r\ is determined. 

In general, there are fewer results available for the directed 
percolation than there are for the undirected percolation. For 
example, the universal exponents Q may be different. Hence, 
in Section IV we investigate the percolation threshold of the 
directed model first by examining the CDF of the cluster size, 
and then observe that the universal exponents may well be the 
same for both cases. To this end, we have developed a fast 
Monte Carlo simulation tool that provides us with samples 
from the cluster size distribution. The pseudo code of the core 
routine is given in Algorithm [T] 

The algorithm partitions the space-time into unit boxes as 
illustrated in Fig. [5] The dimensions are chosen in such a 
way that a node appearing to the shaded box in the center can 
interact during its lifetime only with nodes for which the initial 
position was in one of the 27 surrounding boxes including the 
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Figure 5: Partitioning of the 3D space-time. 



box in the center. This gives a significant reduction in the 
running time of the algorithm, where for each arriving node 
one needs to check the contact times with all other nodes. 
As we may be dealing with infinite clusters, an additional 
condition for termination is needed. To this end, we have 
defined the maximum cluster size s max we are interested in. 

Algorithm 1 Pseudo code for cluster size samples. 

iterate(cylinder c, time t) 

Populate the 27 boxes around c (if not yet populated) 
Store the boxes to hash H for a fast lookup 
for each cylinder c within the neighborhood (27 boxes) do 
if c' overlaps with c then 

r = {z | z 6 overlap} {r is the contact time-interval} 
if max t > t then 

t' — min{z 6 t | z > (} 

if c' is not marked at time t' then 

Schedule event (c, c' , t') at time t c for the pair (c, c') 
end if 
end if 
end if 
end for 

Main routine: 

Scale the dimensions so that a cylinder fits inside a unit box 

Partition the 3D space to unit boxes: (x, y, z) £ Bi.j.k i < 

x < i + 1, j < y < j + 1, k < z < k + 1. 

Populate Bc^o.o an d store it to hash H for fast lookups 

Place the source node Co at (0.5, 0.5, 0) in box Bo,o,o 

Set current time t = and cluster size s = 1 

Set the message acquisition time of node Co to t 

iterate(co, t) 

while scheduled events do 

Pop the next event (c, c' ,t') from a priority queue, set t — t' 
if c' not marked and s < s max then 
s = s + 1 

Set the message acquisition time of node c' to t 

iterate(c' , t) 
end if 
end while 
return s 



C. Asymptotic behavior 

Before proceeding with the numerical results according to 
the percolation theory, let us first discuss the asymptotic behav- 
ior when the distance I is much longer than the transmission 
range d, £ ^> d. In this case, the situation is similar to the 
one analyzed in |2j, pT| : contacts with other nodes are point 
contacts, which occur at a constant rate of 

\ a 8 _ , 
A = — nd, 

7T 



per unit distance on the (x, y)-plane. Near the criticality 
threshold, only a small fraction of nodes carry the message. 
Assuming a node meets at most one node with the message, 
the contact time is uniformly distributed on interval (0,^) 
along the path. In order for the system to avoid extinction, 
each node acquiring the message should pass it further to at 
least one other node on average. That is, 
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s — fidds > 1, 



giving the fluid bound (fb), 



*>*f>(7)^, 



(5) 



which is valid when 7 tends to infinity. At this bound, X-£ = 2, 
i.e., a node meets on average only two nodes before departing. 
Conversely, with aid of |2|, we find a scaling law for the tail 

of Vc{l), 

rj c (^) oc I/7, for 7 > 1. (6) 

It is straightforward to generalize the above to non-constant 
transitions, i.e., to the case where the movement £ is an 
i.i.d. random variable, £ ~ X, with PDF f(x). The point 
contact assumption means that the transitions must be (almost 
surely) such that £ ^> d. Moreover, in a microscopic scale 
the movement is assumed to be isotropic and direct, while in 
a macroscopic scale the paths can make, e.g., smooth turns 



1 21]. The criticality condition is obtained as in above: each 



node acquiring a message should, on average, pass it to one 
new node. According to the theory of renewal processes (see, 
e.g., [22]), the PDF for transition length £ on condition that a 
node obtains a message is 

g(x) = xf(x)/E[X]. 

(cf. the hitchhiker's paradox). The remaining transition length 
R after the acquisition on condition that £ — x is uniformly 
distributed on (0,x) and E[R] = x/2. The number of new 
contacts A in R is 

4 

E[A I £ = x] = -ndx. 

7T 

Hence, the mean number of contacts after the acquisition is 

/>oo 

E[A] = / g(x) E[A | x] dx, 



(7) 



and requiring ~E[A] > 1 then yields, 

ttE[X] 
4E[X 2 ]' 



<,d > 



The second moment in the denominator underlines the fact 
that varying transition length improves the transport capacity 
in DTN communication. 
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Figure 6: Numerical results with stationary nodes, 7 = 0. Note the linear scale on the left figure and logarithmic scale on the right figure. 
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Figure 7: Closer look at the numerical results with stationary nodes, 7 = 0. 



IV. Numerical results 

We have carried out a large number of Monte Carlo simula- 
tions with the aforementioned simulation tool in order to de- 
termine numerical values for rj c (^). The sample size for each 
(77, 7)-pair is 60000 clusters. The parameter s max , defining the 
maximum cluster size, was chosen to be 2 20 + 1 = 104 8 5 77, 
i.e., about 1 million. With these values, the running time is still 
reasonably short with a standard PC. Such clusters are rather 
large creatures, e.g., there are more than 10 12 node pairs, and 
a naive implementation would choke on them. 

A. Stationary nodes 

First we assume no mobility, i.e., 7 = 0, and vary the 
reduced number density, 77 = 0.365, 0.383, which cor- 
responds to the mean node degree of v — 1.46, 1.53. 
If v sa 4.51 or higher, Gilbert's disc model percolates and 
a gigantic component emerges (at any given time instant) 
enabling multi-hop communication (Tl. However, in our case 
the node density is clearly below that and one has to resort to 
DTN-style communication. This is illustrated in Fig. [8] where 
v = 1.52 in the left figure and v = 4.51 in the right figure. 

Recall that the percolation threshold r\ c corresponds to the 
smallest node density at which a gigantic component emerges, 



Therefore, if CDF of the cluster size S remains asymptotically 
below 1, lim s ._ i . 00 P{S < s} = p and p < 1, then the given 
system is above the percolation threshold. The quantity 1 — p, 
referred to as the strength, is the probability that a cluster 
starting from a random node is infinite (TTJ. 

The simulation results are illustrated in Fig. [6] The left 
figure depicts CDF with a linear scale on the x-axis, and the 
right figure with a logarithmic scale. The logarithmic scale 
gives more insight to the behavior. The CDF for 77 = 0.365 
converges to 1 at s « 10 4 , i.e., in practice every realization 
goes to extinction and no message will be delivered for more 
than to about 10 4 nodes. However, the CDF for 77 = 0.379 
seems to converge to a finite value less than 1, suggesting that 
a small fraction of messages may survive forever. Increasing 
77 further to 0.383 improves the chances of hitting a gigantic 
cluster to about 20%, i.e., with rj = 0.383 the strength of the 
percolation is about 0.2. 

Conservatively, we observe that 77 = 0.383 is above the 
percolation threshold, while 77 = 0.375 is most likely below 
it, 0.375 < rj c < 0.383. Comparing the directed percolation to 
the undirected one, we see that the direction constraint causes 
77c to increase from 0.3312(1) to about 0.38, i.e., about 15% 
increment. This is in agreement with the observations made in 



Vc 



inf{?7 I P{S = 00} > 0}. 



1 14], where an actual content dissemination system (obeying 
the causality) was also evaluated by means of simulations. 
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Figure 9: Numerical results with mobile nodes, 7 = 1. 
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Figure 10: Closer look at the results with mobile nodes, 7 = 1. 
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Figure 8: Snapshot of two networks. 



1) Universal exponents: As mentioned, it is not clear 
whether the universal exponents Q are the same for the 
directed model and the normal undirected model. Fig. [7] (left) 
shows CDF of the cluster size distribution in the interesting re- 
gion near the percolation threshold. The graph clearly suggests 
that 0.378 < 77 c (0) < 0.380 and our estimate is ?y c (0) 0.379. 
In Fig. [7] (right), the universal exponents according to Q are 
used. The curve corresponding to 77 = 0.379 has a constant 
tail, which suggests that the universal exponents are the same 
for the undirected and directed percolation models, and we 
assume so in the rest of the paper, (cf. also the Appendix). 



B. Mobile nodes 

Next we set parameter 7 to 1, i.e., each node moves a 
distance equal to the transmission range before departing. The 
numerical results are given in Fig. [9] where the left figure is 
again in linear scale and the right in logarithmic. We observe 
that the critical r? c (l) s» 0.282, as the curve with 77 = 0.279 
clearly converges to 1 and the curve with 77 = 0.286 stabilizes 



at about 0.8. Fig. 10 takes a closer look at the interesting 
regime. Left figure depicts the CDF, which suggests that 
0.282 < T) c (l) < 0.283. In the right figure, we again use 
the universal exponents and the curve with 77 = 0.2825 indeed 
appears to have a constant tail. However, a closer inspection 
of the curve reveals that in fact it is slightly increasing already 
at s a = 500, suggesting ?7 C < 0.2825. As three digits is more 
than enough for our purposes, we conclude that r] c (l) rj 0.282. 

Recall that v — 4 77 according to d2l, i.e., the minimum 
feasible mean node degree v c (j) is obtained by multiplying 
the corresponding 770(7) by 4. According to our numerical 
results, ?7 C (1) ~ 0.75 • ?7c(0), i.e., 25% smaller. This means 
that both the node density n and the mean node degree v can 
also be 25% smaller, i.e., even a small mobility equal to the 
transmission range improves the message forwarding capacity 
considerably (in space-time). 

The critical reduced number density 7^(7) for several other 
values of 7 is determined in the Appendix. Fig. [TT] depicts the 
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Figure 11: Critical mean node degree v c {~i) as a function of the tilt 
ratio 7 characterizing the relative mobility in the model. The fluid 
bound curve is according to |5} and valid at the limit 7 — > 00. The 
other curve is due to percolation theory, f c (j) = 4 r) c (7). 

behavior of the critical mean node degree ^(7) as a function 
of 7. The fluid bound curve is according to |5]l and valid at 
the limit 7 — > 00. The other curve corresponds to the critical 
percolation threshold, v c {f) = 4 770(7), the values of which we 
have obtained numerically in this paper. We observe that the 
highest gain from an additional movement of A£ is obtained 
when 7 = 0, as expected. 

The two curves differ initially, where |5} is not valid. 
However, already at 7 w 4 the two curves are almost matching. 
Consequently, for a small movement one needs to analyze the 
situation according to the percolation, but as soon as, say, 
7 > 4, the fluid flow analysis already gives accurate results. 

V. Conclusions 

DTN is designed to operate in settings where connectivity is 
intermittent at best. In this paper, we have analyzed necessary 
conditions for DTN type communication. To this end, we 
assumed an elementary mobility model where nodes join the 
network for a constant time t and then depart. During the 
sojourn time, they move a distance £ to a random direction. 
We assume a large network, where the ability to sustain a 
message indefinitely means that the message will eventually 
also reach its final destination(s). The stochastic model of DTN 
is amenable to analysis by means of the percolation theory. 
In particular, we studied the so-called directed continuum 
percolation in space-time, where the objects are cylinders with 
height t and the diameter equal to the transmission range d 
sheared (tilted) around the time axis by an amount correspond- 
ing to the movement £. We showed that the critical percolation 
threshold in terms of the critical reduced number density r\ c 
depends only on the tilt ratio 7 = Ijd, r\ c = 770(7), where 
I denotes the distance nodes move and d is the transmission 
range. Numerical values for 770(7) were obtained by means of 
Monte Carlo simulations, where universal exponents of the 
undirected model were utilized. Some evidence to support 
the hypothesis that the universal exponents are the same for 
the directed model was also given. Moreover, the asymptotic 



behavior of 770(7) when 7 tends to infinity was derived. In 
terms of the mean node degree v, our main result states 
that v > 4770(7) m order for a large DTN network to be 
operational. A general observation is that increased mobility 
improves the transport capacity of DTN and allows a lower 
node density, as expected. 

The elementary mobility model, as well as, the symmetric 
Gilbert's disc communication model, represent, in some sense, 
a strict lower bound for more realistic scenarios. That is, if, 
e.g., the sojourn time t or the distance i were random variables, 
or, if directional antennas were used, then the shapes and 
the orientation of the objects would vary more, increasing 
"disorder" and yielding a lower critical percolation threshold. 
Such studies are left as a future work. 
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7 vein) 


7 ndn) 


0.0 0.3788 
0.05 0.3735 
0.1 0.368 
0.2 0.357 
0.5 0.3262 


1 0.2825 

2 0.2189 
4 0.1459 
8 0.0836 

16 0.04325 



Table I: Critical reduced number density r? c (7) based on the Monte 
Carlo simulations, where ri c (8) ~ 2?7 C (16) in agreement with {6|. 



Appendix 

An estimate for the critical reduced number density r/ c (j) is 
determined for 7=0.05, 0.1, 0.2, 0.5, 2, 4, 8, 16 in Figs. [12}- 
13 In each case, parameter 77 obtains values from both sides 
of the corresponding critical percolation threshold. Table [I] 
summarizes the results. In Fig. 14 the maximum cluster size 
is very large, s max = 10 8 , suggesting 7i c (0) = 0.3788(1) and 
also supporting the hypothesis on the universal exponents. 
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Figure 14: Results with 7 = and s m ax = 10 s , i.e., 100 million. 



